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£f) ' Abstract. - We present a new framework for modeling the statistical behavior of both fully developed 

turbulence and short-term dynamics of financial markets based on the generalized non-extensive thermo- 
statistics formalism. We also show that intermittency - strong bursts in the energy dissipation or clusters of 
high price volatility - and non-extensivity - anomalous scaling of usually extensive properties like entropy - 
are naturally linked by a single parameter q, from the non-extensive thermostatistics. 

Ctf ; 

Scaling invariance plays a fundamental role in many natural phenomena and frequently emerges 
■ from some sort of underlying cascade process. A classical example is fully developed homogeneous 

isotropic three-dimensional turbulence, which is characterized by a cascade of kinetic energy from 
large forcing scales to smaller and smaller ones through a hierarchy of eddies. At the end of 
the cascade, the energy dissipates by viscosity, turning into heat. Recently, some authors [1-3] 
have studied the phenomenological relationship among financial market dynamics, scaling behavior 
and hydrodynamic turbulence. Particularly, Ghashghaie et al. [2] conjectured the existence of 
a temporal information cascade similar to the spatial energy cascade found in fully developed 
turbulence. 

Traditionally, the properties of turbulent flows are studied from the statistics of velocity differences 
v r (x) = v(x) — v(x + r) at different scales r. As with other physical systems that depend on the 
dynamical evolution of a large number of nonlinear ly coupled subsystems, the energy cascade in 
turbulence generates a spatial scaling behavior - power-law behavior with r - of the moments (v™) 
of the probability distribution function (PDF) of v r (the angle brackets () denote the mean value 
of the enclosed quantity). For large values of the Reynolds number, which measures the ratio of 
nonlinear inertial forces to the linear dissipative forces within the fluid, there is a wide separation 
between the scale of energy input (integral scale L) and the viscous dissipation scale (Kolmogorov 
scale rj). Though at large scales (~ L) the PDFs are normally distributed, far from the integral scale 
they are strongly non-Gaussian and display wings fatter than expected for a normal process. This 
is the striking signature of the intermittency phenomenon. After publication of the Kolmogorov 
K62 refined similarity hypotheses [4], the problem of small scale intermittency became one of the 
central questions on isotropic turbulence. Over the past years several papers [5-12] have discussed 
intermittency and the so-called 'PDF problem'. Similar attempts [1,2,13] have been made to explain 
the same peculiar shape observed in the PDF of price changes z T = z(t) — z(t + r) at small time 
intervals. 
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Based on the scaling properties of multifractals, a generalization of Boltzmann-Gibbs thermo- 
statistics has been proposed [14-17] through the introduction of a family of non-extensive entropy 
functionals Snip) with a single parameter q. These functionals reduce to the classical, extensive 
Boltzmann-Gibbs form as q — > 1. Optimizing S q \p] subject to appropriate constraints [16], we 
obtain the distribution 

p q (x) = [l-P(l-q)x 2 ] l /^/Z q . (1) 
The normalization factor, for 1 < q < 3, is given by 

V 3 r((3 - q)/2(q - 1)) 
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In the limit of q — > 1, we recover the Gaussian distribution. 

The above distribution, we claim, provides a simple and accurate model for handling the PDF 
problem. To show this, we stay in the context of fully developed turbulence (x = v r ). From 
equation (|I|), we can easily obtain the second moment 

and the flatness coefficient 

K _{4)_ _ 3(5-3g) 
r "(^ 2 ) 2 " (7-5,) • [6) 
We remark that the flatness coefficient, which is directly related to the occurrence of intermittency, 
is solely determined by the parameter q. 

At this point, if we assume [2,5,6] a scaling of the moments (v™) of v r as r^ n , the variation with r 
of the PDF of the velocity differences and of its related moments can be completely determined. 
Particularly, we can obtain the functional forms of the flatness coefficient and the parameter q, 
respectively 

K r = K L (r/L) a (4) 

and 

= lh-7K L {r/L) a = 15 - 7K v {r/r)) a 
q 9-5K L (r/L) a 9-5K v (r/r)) a ' [ ' 

where K v is given by equation (||), a = — 2C2 and Kl = 3, the expected value for a Gaussian 
process. The correspondent expression for j3 can be derived similarly from equation (0). For a PDF 
normalized to its standard deviation (that is, with unit variance), we have (3 = 1/(5 — 3q). 

We note that in the limit of infinite Reynolds numbers, as r — > 0, K r diverges (since a < 0) while 
q tends to a finite limit, q < 7/5. This bound coincides with the one obtained by Boghossian 
[18] through a ^-generalization of Navier-Stokes equations. Moreover, this limit implies that the 
second moment of distribution (||) will always remain finite, which is empirically expected from the 
phenomena here analyzed. 

In order to account for the well known asymmetry of the velocity distributions, we may also consider 
P = P p , for v r > 0, and P = p m , for v r < 0. In this case, we have 

K)=MPrn,P P )s^ L -\^f 2 j , (6) 

r (o) r ( S - 2) 



where s = l/(g — 1) and 



( n+ l)/2 , s n ( n+ i)/2 
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Accordingly, the value of parameters q, (3 m and (3 p , at a given scale r, are now determined from the 
second moment, the skewness (S r = (vf) / (i^) 3 / 2 ) and the flatness coefficients. 

We checked our model with turbulence statistics data taken from reference [2] , provided by Chabaud 
et al. [9]. For this, we compared the experimental data with the theoretical predictions, for an 
skewness factor at the integral scale of Sl = —0.4, L = 1cm [9] and a = —0.10 [5,19]. The results 
are displayed in Fig. la. A good agreement is observed through spatial scales spanning two orders 
of magnitude - from the neighborhood of the integral scale down to close to the Kolmogorov scale 
— , and for a range of up to 15 standard deviations, including the rare fluctuations in the tails of 
the distributions. Note that the solid lines in Fig. la have not been adjusted to the data through 
a free parameter, as in other models [9,6,7]. In the present case, the parameters q, (3 m and (3 P , that 
control the shape of the PDF in each scale, are uniquely determined from the scaling of (v^), S r 
and K r , obtained from equation (||). 

The same approach adopted in turbulence can be straightforwardly applied (with x = z T and 
z T scaling as r^ n ) to model the statistics of price differences in financial markets, as far as the 
parameter a, Sl and the integral time scale tl ~ time span for which a convergence to a Gaussian 
process is found - are available. We tested our model with price changes data taken from reference 
[2], provided by Olsen & Associates. The results are displayed in Fig. lb, for a = —0.16, Sl = —0.4 
and tl — 2.3 days. Again, we observe that the proposed model reproduces with good accuracy the 
statistics of price differences over all temporal scales. 

Non-extensivity, a matter of speculation in some areas [20], is an essential feature of the generalized 
thermostatistics. If we suppose a scenario of a cascade of bifurcations with m levels, and scale the 
generalized expectation value of an observable O q (the kinetic energy ^u 2 of velocity differences, 
for example), averaged over a volume of size V = rf and normalized by Boltzmann constant, we 
have at the first level [21] 

O g (2V) = 20 q (V) + 2(1 - q Q )O q (V)S q (V) (8) 

and at level m 

O q {2 m V) = 20 q {2 m - 1 V) + 2(1 - q m ^)O q {2 m - l V)S q {2 m - l V) 

^2 m O g (y) + (l-g m _ 1 )2 m+c O (? (y) , (9) 

where the higher order terms in (1 — q) have been neglected, with C being a constant to be 
determined later. Cascade processes are also described in terms of fractal or multifractals models 
[7,22-25]. Within these frameworks, in high Reynolds number turbulence, the energy dissipation 
is not uniformly distributed within the fluid but rather concentrated on subsets of non-integer 
fractal Dp dimension. This picture leads to a scaling behavior with dimensionality not equal to 
the dimension D of the embedding space. In this case, if we consider the cascade of bifurcations 
described above, we find 

O q (2 m V) = 2 mDp /D O q (V) , (10) 



with D = 3. It follows immediately from equations (^) and ([!(]) that 

~log(2- c + 1 - 
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At the top of the cascade, we have q m -i = 1 an d Dp = D. C is determined from the value of Dp 
at the bottom of the cascade. Writing 2 m ~ 1 r/ = r and g m _i = q(r), we get 

^ rt log((2- c + l- g (r))2 c + 1 r/>7) 

D F ~D — , (12) 

log(2 r/r/J 

where g(r) is given by equation @. For C = —1, at the bottom of the cascade (r = 77), using the 
values of a and L previously specified, and 77 = 0.022mm [9], we get Dp ~ 2.37. This value is 
in good agreement with the fractal dimension of interfaces in turbulent flows (Dp = 2.35 ± 0.05), 
measured in different experimental contexts [26-29]. 



Equation (12), through the variation of parameter q, offers a quantitative picture of the transition 
from small-scale intermittent, non-extensive, fractal behavior to large-scale Gaussian, extensive 
homogeneity. This equation can also be exploited to estimate the structure functions exponents. 
Under the assumptions of the multifractal model [30], Q n is given by Q n = inf/Jn h+D— Dp(h)], with 
(3 = 1, assuming that the turbulent flow possesses a range of scaling exponents I = (h m i n , h max ). 
The existence of a corresponding range of dissipation scales if extending from Tj m i n *~ 
to 

Vmax ~ L R V(i+^ ma:c ) ; where R is the Reynolds number at the integral scale, allows us to rewrite 

C n = inf [nhtf)+D-D F tf)] , (13) 



where Dp(rj') is given by equation (|12|), with r = 77'. The resulting exponents £ n are in good 
agreement with experimental values [19,5], as shown in Table 1, for n = 2, . . . , 8. 

The above picture may also be applied to the information cascade, with D = 1. One main qual- 
itative difference between the two processes is that, since there is nothing equivalent to viscous 
damping in the dynamics of speculative markets, the information cascade depth is only limited by 
the minimum time necessary to perform a trading transaction (roughly 1 min [1]). On the other 
hand, as in turbulence, the scaling exponents depend on the order of the moments in a nonlinear 
way. Figure 2 displays the comparison of the experimentally measured scaling exponents £ n [2] and 
the theoretical prediction. We found a good coincidence between our results and the experimental 
data, assuming an equivalent Reynolds number R of approximately 9000 and £2 = 1-06 [3]. 

Summarizing, we described a simple and accurate framework for modeling the statistical behavior 
of both fully developed turbulence and short-term dynamics of financial markets based on the 
formalism of the generalized non-extensive thermostatistics. Within this framework, we have shown 
that intermittency and non-extensivity are naturally linked by parameter q, which represents an 
objective measure of small-scale intermittent, fractal behavior in turbulent cascades. 
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Captions 

Table 1 Comparison of experimentally measured scaling exponents £ n and the theoretical predic- 
tion. 

Figure 1(a) Data points: standardized probability distributions p q {v r ) of velocity differences 
v r (x) = v(x) — v(x + r) for spatial scales, from top to bottom, r/77 = 3.3, 18.5, 138 and 325 (data 
taken from ref. [2], provided by Chabaud et al. [9]); Solid lines: proposed PDF model, equation 
(|l|) with j3 = /3 P , for v r > 0, and /3 = /3 m , for v r < (for better visibility the curves have been 
vertically shifted with respect to each other). 

Figure 1(b) Data points: standardized probability distributions p q (z T ) of price differences z T {t) = 
z(t) — z(t + t) for temporal scales, from top to bottom, r = 640 s, 5120 s, 40960 s and 163840 s 
(data taken from ref. [2], provided by Olsen & Associates); Solid lines: proposed PDF model, 
equation (||) with f3 = (3 p , for z T > 0, and (3 = (3 m , for z T < (for better visibility the curves have 
been vertically shifted with respect to each other). 

Figure 2 Comparison of theoretical and experimentally measured scaling exponents £ n . 
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